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ABSTRACT 

We present the broad-band noise structure of selected Anomalous X-Ray Pulsars 
(AXPs) and Soft Gamma Repeaters (SGRs) in the 2-60 keV energy band. We have an- 
alyzed Rossi X-ray Timing Explorer Proportional Counter Array archival light curves 
for four AXPs and one SGR. We detect that the persistent emission of these sources 
show band limited noise at low frequencies in the range 0.005-0.05 Hz varying from 
2.5% to 70% integrated rms in times of prolonged quiescence and following outbursts. 
We discovered band-limited red noise in IE 2259+586 only for ~2 years after its major 
2002 outburst. The system shows no broad-band noise otherwise. Although this rise 
in noise in IE 2259+586 occurred following an outburst which included a rotational 
glitch, the other glitching AXPs showed no obvious change in broad-band noise, thus it 
does not seem that this noise is correlated with glitches. The only source that showed 
significant variation in broad-band noise was IE 1048.1-5937, where the noise gradu- 
ally rose for 1.95 years at a rate of ~3.6% per year. For this source the increases in 
broad-band noise was not correlated with the large increases in persistent and pulsed 
flux, or its two short SGR-like bursts. This rise in noise did commence after a long 
burst, however given the sparsity of this event, and the possibility that similar bursts 
went unnoticed the trigger for the rise is noise in IE 1048.1-5937 is not as clear as for 
IE 2259+586. The other three sources indicate a persistent band-limited noise at low 
levels in comparison. 
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1 INTRODUCTION 

Anomalous X-Ray Pulsars (AXPs) and Soft Gamma Re- 
peaters (SGRs) are part of a new class of neutron stars 
(NSs), whose emission mechanisms commute with neither 
accreting X-Ray Pulsars nor radio pulsars (PSRs). Their 
periods are clustered between 2 - 12 s and they are charac- 
terized by their rapid spin down (~ 10~ 13 — 10 -11 s s _1 ) and 
high quiescent luminosities (~ 10 33 — 5 x 10 35 erg s _1 ) that 
can not be explained by spin down mechanisms. They also 
show wide range of variability; including flux and pulse pro- 
file variations (Dib, Kaspi & Gavriil 2007a; Tam et al. 2008), 
glitches (Dib, Kaspi & Gavriil 2007b), bursts (Gavriil, Kaspi 
& Woods 2002) and strong timing noise (Baykal et al. 2000; 
Kaspi, Gavriil & Chakrabarty et al. 2001). Three SGRs are 
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known to exhibit giant bursts (SGR 0526-66, SGR 1900+14, 
SGR 1806-20), which in two of these objects rapid quasi- 
periodic oscillations (QPOs) are detected after the onset of 
their flares (SGR 1900+14; Strohmayer & Watts 2005, SGR 
1806-20; Israel et al. 2005). 

SGRs and AXPs are widely accepted as being mag- 
netars, neutron stars (NSs) with magnetic fields strengths 
greater than the quantum critical level (10 14-15 Gauss), 
whose main source of energy is the dissipation of these fields 
(Duncan & Thompson 1992; Thompson & Duncan 1995). 
The magnetar model for these objects explains the SGR 
bursts and the large flares, and it predicted the AXP bursts 
which occurred in IE 1048.1-5937 (Gavriil, Kaspi & Woods 
2002) and IE 2259+586 (Gavriil & Kaspi 2002). Thompson 
& Duncan (1995) argue that the bursts are due to crustal 
cracking as a result of magnetic diffusion, and that the large 
flares are due to major reconfigurations of the global mag- 
netic field (Thompson & Duncan 1995). Magnetar strength 
magnetic fields can also account for the rapid spin down of 
these objects and the QPOs observed during large flares (e.g, 
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Israel et al. 2005; Strohmayer & Watts 2005). The existence 
of a fossil disk, and accretion from it, has also been pro- 
posed to explain the persistent properties of these sources 
(van Paradijs et al. 1995; Chatterjee et al. 2000; Alpar 2001). 
However, this model has difficulty explaining the giant flares, 
repeat bursts and the pulsed optical emission that exists in 
this class of objects. 

In this paper, we investigated the properties of the 
broad-band noise in selected AXPs and SGRs, and present 
plausible origins for it within the context of the above two 
models. For example, if AXPs and SGRs have fall-back disks 
that they accrete from, as has been suggested by Ertan et 
al. (2007), then we can compare this noise to that of other 
accreting sources. For instance, the broad-band noise in X- 
ray binaries spans a very broad range and is due to variable 
X-ray flux from the accretion disk. In corona hosting X-ray 
binaries, a contribution of the noise in the high frequency 
range is due to photon delays in Compton up scattered ra- 
diation (van der Klis 2006 and references therein). This is 
important given that AXPs and SGRs show hard power law 
like tails in their energy spectra that are also most likely 
due to photons being up-scattered in energy. In the context 
of the magnetar model, these hard-tails in AXPs and SGRs 
have been successfully modeled by resonant cyclotron scat- 
tering (Lyutikov & Gavriil 2006; Rea et al. 2008). We argue 
that resonant cyclotron scattering can also introduce broad 
band noise in the frequency range observed for these sources. 
Thus, in order to discriminate between the two models it 
becomes of great interest to study the broad-band noise in 
AXPs and SGRs. 



2 DATA AND OBSERVATIONS 

To calculate the band limited noise of these objects we used 
the 2-60 keV archived public Rossi X-ray Timing Explorer 
(RXTE) light curves. The RXTE archive has a total of 10 
sources that are either AXPs or SGRs. We used only ap- 
propriate sets of data, eliminating sources depending on 
epochs, exposures, burst dates and objects in the field of 
view (FOV). 

We recovered 1 known SGR and 4 known AXP ob- 
servations with fields relatively free of other sources. The 
exceptions are IE 1048.1-5937 with Eta Carinea and SGR 
1806-20 with XTE J1810-197 in their FOV, which both do 
not effect the noise levels in our analysis. The variability 
of Eta Carinae occurs in a very different time scale (~85 
days; Corcoran et al. 1997) with respect to 256 seconds 
used in this study. XTE J1810-197 is a transient source 
which was quiescent during the observations in our analysis 
(Ibrahim et al. 2004). Other fields include strong and vari- 
able sources. The data were obtained by the Proportional 
Counter Array (PCA; Jahoda et al. 1996) instrument aboard 
RXTE. We obtained 125 ms resolution background sub- 
tracted light curves from archived standard products (Std- 
Prods) of RXTE. These light curves are screened through 
standard procedures; namely source elevation greater than 
10°, pointing offset smaller than 0.02°, also PCU turn- 
off times are subtracted through defining good time inter- 
vals. Backgrounds are generated by pcabackest command of 
FTOOLS with the use of background models of the relevant 
epoch. These net lightcurves also include the barycentrically 



Table 1. Observation log of the light curves analyzed in this 
work. 



Object Obs Id Ave. Count Rate Burstt 



c/s 


IE 1048.1-5937 


40083- 


-08 


3.889 


no 


IE 1048.1-5937 


50082- 


-04 


3.771 


no 


IE 1048.1-5937 


60069- 


-03 


3.302 


yes 


IE 1048.1-5937 


70094- 


-02 


7.030 


no 


IE 1048.1-5937 


80098- 


-02 


5.667 


no 


IE 1048.1-5937 


90076- 


-02 


3.580 


yes 


IE 1048.1-5937 


91070- 


-02 


3.255 


no 


IE 1048.1-5937 


92006- 


-02 


3.457 


no 


IE 2259+586 


20145- 


-01 


1.799 


no 


IE 2259+586 


20146- 


-01 


2.269 


no 


IE 2259+586 


40082- 


-01 


1.752 


no 


IE 2259+586 


50082- 


-01 


1.167 


no 


IE 2259+586 


60069- 


-01 


0.4600 


no 


IE 2259+586 


70094- 


-01 


2.695 


yes 


IE 2259+586 


80098- 


-01 


1.765 


no 


IE 2259+586 


90076- 


-01 


1.673 


no 


RXS J1708-40 


40083- 


-11 


13.02 


no 


RXS J1708-40 


50082- 


-07 


12.10 


no 


RXS J1708-40 


60069- 


-07 


12.37 


no 


RXS J1708-40 


60412- 


-01 


12.99 


no 


RXS J1708-40 


70094- 


-01 


13.58 


no 


RXS J1708-40 


80098- 


-04 


13.38 


no 


RXS J1708-40 


90076- 


-04 


13.38 


no 


IE 1841-045 


40083- 


-11 


9.521 


no 


IE 1841-045 


50082- 


-05 


8.835 


no 


IE 1841-045 


70094- 


-03 


9.937 


no 


IE 1841-045 


90076- 


-03 


9.013 


no 


IE 1841-045 


91070- 


-03 


9.006 


no 


SGR 1806-20 


50142- 


-01 


7.626 


yes 


SGR 1806-20 


50142- 


-03 


7.607 


yes 


SGR 1806-20 


60121- 


-01 


7.053 


yes 


SGR 1806-20 


70136- 


-02 


10.69 


yes 


SGR 1806-20 


90073- 


-02 


10.93 


yes 


SGR 1806-20 


90074- 


-02 


12.70 


yes 


SGR 1806-20 


91062- 


-02 


8.721 


yes 



t Observation dates have been correlated with the GCN 
(http://gcn.gsfc.nasa.gov/gcn3_archive.html) burst data 
and relevant papers (see captions of Fig. 5). 



corrected time columns, which were used in our analyses^. 
Only when necessary, particularly to check background ef- 
fects, non-background subtracted light curves were created 
from original data using SEEXTRCT v4.2. The manipula- 
tion of the data was made with the FTOOLS v5.21 software. 

In order to measure the broad-band noise, we derived 
and averaged several power spectra for each source. The 
power spectral densities (PSD) expressed were calculated 
in terms of the fractional rms amplitude squared following 
Miyamoto et al. (1991), and the expected white noise levels 
were subtracted hence leaving us with the rms fractional 
variability of the time series in units of (rms/mean) 2 /Hz. 
We used in total 658 light curves for our analysis (151 for 
IE 2259+586, 136 for IE 1048.1-5937, 125 for SGR 1806- 
20, 82 for IE 1841-045 and 153 for RXS J1708-40), each of 



1 All of the procedures involved in or- 
der to produce StdProds are explained in 
http://heasarc.gsfc.nasa.gov/docs/xtc/rccipcs/stdprod_guide.html 
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Figure 1. The averaged PSD in time showing the broad- band 
noise of AXP IE 1048.1-5937. MJD 52211 and 52227 arc burst 
epochs (Gavriil & Kaspi 2004). The time line increases top to 
bottom first, then from left to right. 



which was segmented into 10-50 individual power spectra 
to construct average PSD in a given time line. 



3 ANALYSIS AND RESULTS 

In order to increase the signal to noise of the PSD to look 
for broad-band noise, the averaged light curves are grouped 
in three month intervals yielding 20 ks-50 ks long exposure 
times for each source. The burst intervals noted in Table 1 
are excluded from our data samples to avoid additional red 
noise because of the light curve variability during bursts. The 
remaining burst like features ( especially for the bursting SGR 
1806-20) and spurious events due to instrument errors are 
screened by hand from the light curves. We averaged several 
256 s long PSD for each source. 

Figures 1 through 4 show the representative av- 
eraged PSD fitted with a broken power law model; 
7 [l. + (x/u) 2 ~^ " , in which v is the frequency, 7 is the 
(rms/mean) 2 /Hz value of the flat part, x is the cutoff fre- 
quency and a is the power law index. Each spectra in the 
analyses is averaged over three-month intervals and are sep- 
arated by 5-900 days. In all sources, we find red noise with 
a break in the averaged PSD which may be considered band 
limited noise (BLN). 

For all figures, the reduced chi squares of the fits are 
between 0.6-1.2, the range of observed break frequencies are 
between 0.05-0.7 Hz and the power law indices are between (- 
)0. 5-375. Table 2 shows the list of power law indices derived 
from the fits to PSD of the sources in time. An irregular 
variation of the indices in time is apparent. The ranges in 
Table 2 denote 1 a confidence level error ranges from fits. 
It is very hard to determine the value of these indices, since 
the decline in the power spectrum is near the edge of the 
sampling window and errors also increase considerably. We 
would like to make a note that for 2 a errors on the power law 
index, the canonical value of -2 is in the acceptable range. In 
addition, the fitted PSD are integrated over 5 x 10~ 3 - 5 x 
10~ 2 Hz interval to get the integrated fractional rms values. 



Time evolution of these values for IE 1048.1-5937 and IE 
2259+586 are plotted in Figure 5. 

We found that for IE 1048.1-5937 and IE 2259+586 
(two AXPs), their burst (maybe also flare) epoch correlate 
with their broad-band noise levels. The steady ~ 15% rms 
noise of IE 1048.1-5937 showed a slight drop, during which 
two of its SGR-like short bursts and a flare appeared, to 
levels around ~ 10% rms in 16 months. This low noise re- 
covered throughout the next 5 months until it reached its 
former noise level (see Figures 1 and 5). The red noise level of 
the source rose after a long burst and later a glitch and large 
scale flux increase (Dib et al. 2007; Tam et al. 2008) to levels 
~ 25% rms and stayed at this level (with some variability) 
for about 1.95 years (see Fig. 5, top panel). IE 2259+586 
showed a slightly different behavior. It had no detectable 
broad-band red noise in quiescence and developed broad- 
band red noise around ~ 60% rms level right after its burst 
active state and the on-set of a glitch, stayed consistently 
at this level throughout nearly 2 years, until it diminished 
in about 1.5 months (at MJD 53233) (see Figures 2 and 5 
bottom panel) back to its original level of no broad-band 
noise. 

The other two AXPs IE 1841-045 and RXS J1708-40 
maintain a persistent broad-band noise level compared with 
the two AXPs discussed in the above paragraph. Sources 
showed around ~ 11% rms and ~ 7% rms level noise for IE 
1841-045 and RXS J1708-40, respectively (with slight in- 
crease in time). The total time span of the detected broad- 
band noise for the two different epochs for RXS J 1708-40 
were 2.3 years (51215-52053 MJD) and 2.6 years (52366- 
53325 MJD) . Also for IE 1841-045 the time spans for two 
separate time epochs analized in this work were 2 years 
(51977-51260 MJD) and 3.5 years (53635-52349 MJD). IE 
1841-045 also stands out as the noisier of the two. SGR 1806- 
20 exhibited broad-band noise slightly fluctuating about a 
persistent level of ~ 10%, overall. The time span of our anal- 
ysis in two epoch ranges were 1.4 years (51600-52100 MJD) 
and 1.7 years (53000-53600 MJD). We had difficulty relating 
the burst activity and the PSD for SGR 1806-20 due to the 
existence of excessive number of bursts. 



4 DISCUSSION 

Figure 5, shows the comparative evolution of the band- 
limited red noise in IE 1048.1-5937 and IE 2259+586. The 
ordinate is the integrated rms band-limited noise and the 
abscissa is the time in MJD. The error in time (x error) is 
the total time of the accumulated observation IDs for each 
data point. The y error is the 1 a error in the integrated 
rms values. IE 2259+586 showed the highest level of noise 
and longest duration of change. IE 1048.1-5937 is the only 
source that showed a drop of the noise level assuming level 
at start of MJDs is the steady normal level of noise. It was 
also the most noisy source after IE 2259+586. A significant 
rise of the rms broad-band noise level in IE 1048.1-5937 (like 
IE 2259+586) after the long-burst epoch at 53185 MJD is 
detected for about 1.95 years (may have lasted about 3 years 
judging from telegrams; Rea et al. 2007). 

On the basis of torque noise, IE 2259+586 is very quiet 
(Baykal et al. 2000; Gavriil & Kaspi 2002), but after its 
outburst and a glitch, IE 2259+586 experienced rapid spin 
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Table 2. The power law indices of AXPs and SGRs derived from 
fitted PSD in this work. 





Figure 2. The averaged PSD in time showing the broad- band 
noise in AXP IE 2259+586 with MJD 52443 as the burst epoch 
(Kaspi ct al. 2003). The time line increases top to bottom first, 
then from left to right. 



down over 60 days (Woods et al. 2004) and showed an ex- 
cess in infrared radiation (Hulleman et al. 2001). In this pe- 
riod IE 2259+586 displayed band-limited noise lasting for 
a timescale of 2 years and no band-limited red noise other 
wise. IE 1048.1-5937 shows a highly variable torque noise 
such that, it can not be phase connected for long time in- 
tervals which is consistent with accreting sources (Baykal 
ct al. 2000; Kaspi, Gavriil & Chakrabarty et al. 2001). In 
our work, this source showed band-limited red noise at all 
times. The substantial increase of the broad-band noise of 
IE 1048.1-5937 in Figure 5 (from 15% to about 25-30%) oc- 
cured gradually after the the long burst (> 699 sec) at 53185 
MJD (Gavriil, Kaspi & Woods 2006). Following this, there 
has been a flux increase (outburst-like) from IE 1048.1-5937 
reported along with a spin-up glitch at 54186 MJD (Dib et 
al. 2007; Tam et al. 2008). The decrease of this rise of the 
soft X-ray flux was reported to be at 54265 MJD (Rea et 
al. 2007; Tam et al. 2008). These latest dates are at the 
edge of our curve for IE 1048.1-5937 in Figure 5 and are not 
included in our results. However, it is clear that it follows 
our peak of integrated rms values for the source. We note 
that the rise in broad-band noise started much earlier than 
the flux outburst right after the long burst at 53185 MJD 
( Gavriil, Kaspi & Woods 2006). On the other hand, we cau- 
tion that given the sparsity of this event, and the possibility 
that similar bursts went unnoticed, this start for the rise in 
noise is not as clear as in IE 2259+586. We also note that 
this increase in red noise is not very correlated with a large 
increase in persistent and pulsed flux in the early stages 
until the outburst-like flux increase. Prior to its burst at 
53185 MJD and the broad-band noise increase, there was a 
slight decrease of red noise during which two SGR-like short 
bursts (52211 MJD and 52227 MJD) and two flares (52218 
MJD and 52444 MJD) appeared (Gavriil, Kaspi & Woods 
2002; Gavriil & Kaspi 2004). The two bursts showed char- 
acteristics similar to SGR bursts like marginal pulsed flux 
increase not more than 3cr and short duration (<50 s 1st 
burst and < 2 s 2nd burst). Since we did not detect any cor- 
relation of broad-band noise with the SGR bursts in SGR 



Epoch (MJD) 


Epoch range (d) 


indcx(-) 


index error 


IE 1048.1-5937 








51321 


± 87 


121.97 


± 43.73 


52058 


± 50 


9.99 


± 5.31 


52160 


± 46 


51.64 


± 25.73 


OZZ / 


1 /to 

± 4z 


OZ. ID 


I in ni 

± 1U.U1 


52478 


± 48 


79.31 


± 18.1 


52626 


± 42 


279.25 


± 64.75 


52777 


± 43 


9.24 


± 2.39 


53010 


± 48 


191.55 


± 45.75 


IE 2259+586 








52490 


± 45 


29.07 


± 4.5 


OZDU1 


_|_ pro 


A 77 
4. f 1 


4-917 
± Z.ll 


52757 


± 44 


90.01 


± 38.79 


52863 


± 50 


17.94 


± 3.35 


52975 


± 47 


256.15 


± 65.75 


53087 


± 50 


64.87 


± 27.39 


53192 


± 41 


144.65 


± 34.05 


SGR 1806-20 








51693 


± 76 


0.46 


± 0.08 


51825 


± 46 


0.98 


± 0.2 


52065 


± 44 


16.65 


± 6.41 


53288 


± 43 


1.72 


± 0.19 


53148 


± 52 


108.47 


± 9.74 


0o4yo 


1 A Q 


oZ. 11 


I A go 


IE 1841-045 








51334 


± 74 


18.37 


± 5.29 


51534 


± 63 


6.32 


± 2.06 


51788 


± 60 


24.13 


± 9.54 


51931 


± 46 


78.41 


± 38.69 


52415 


± 66 


30.36 


± 5.78 


52598 


± 69 


30.79 


± 7.52 


52803 


± 77 


124.95 


± 33.76 


52940 


± 40 


45.38 


± 10.03 


53125 


± 52 


39.15 


± 9.16 


53486 


± 46 


211.9 


± 57.5 


53590 


± 44 


13.37 


± 2.47 


1RXS J1708-40 








51317 


± 102 


46.69 


± 16.39 


51530 


± 84 


335.05 


± 179.05 


51845 


± 190 


371.07 


± 282.24 


52044 


± 9 


9.26 


± 4.74 


52409 


± 43 


0.63 


± 0.14 


52538 


± 51 


46.23 


± 11.45 


52681 


± 37 


375.75 


± 172.65 


52794 


± 49 


26.86 


± 8.55 


52906 


± 47 


160.75 


± 30.15 


53009 


± 49 


66 


± 26.16 


53119 


± 35 


24.87 


± 6.61 


53206 


± 47 


9.83 


± 2.04 


53291 


± 34 


53.45 


± 14.31 



1806-20, we do not expect any effects here. The first flare 
starts between the two SGR-like bursts at 52218 MJD, after 
which an epoch of lower rms noise started. The drop of the 
broad-band noise showed a brief increase in the second flare 
(maximum at 52444.4 MJD, see Fig. 5) and the trend of de- 
crease and recovery of noise after the flare is similar to the 
torque variability that was observed by Tam et al. (2008) . 
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The time-averaged PSD of the four AXPs and the SGR 
showed in general very low frequency noise. We detected that 
this noise to be band-limited and in the frequencies higher 
than 0.05 Hz. There were little or no noise at all times be- 
low these frequencies given the sensitivity of RXTE (also 
background) and statistical quality of RXTE light curves. 
The broad-band noise of AXPs and SGRs can be compared 
with the other isolated pulsars or X-ray binaries on the ba- 
sis of their origin (e.g., accretion process, coronal activity, 
rotational variation). Broad-band noise in X-ray binaries is 
a result of aperiodic X-ray flux variability in the form of 
continuum noise and QPOs with frequencies ranging from 
several mHz to more than thousand Hz (Hasinger & van 
der Klis, 1989, van der Klis 2000,2006). The flattening in 
the lower frequencies resemble PSD in disk-fed X-ray bina- 
ries, specifically the low hard state (LS) noise of black hole 
candidates and the very low frequency noise (VLFN) of low- 
magnetic field neutron star LMXBs (van der Klis, 2006 and 
references therein). A flat topped noise arises in black hole 
candidates with large integrated rms noise level (~ 50%) 
in the low frequencies (van der Klis 2006). The VLFN of 
LMXBs show the flat topped noise with a few percent inte- 
grated rms and the continuum power law component shows 
indices between (-)1.5-2, in general a range which holds for 
all the X-ray binaries (van der Klis 2006). The VLFN of the 
X-ray binaries is correlated with M variations associated 
with an accretion disk and plausibly related to unsteady 
nuclear burning (van der Klis 2006). The structure of the 
flattening in the PSD sample of AXPs and SGRs resembles 
the VLFN of accreting X-ray binaries. 

If one assumes that the variable X-ray flux originates 
from an accretion disk with a size between lxlO 9 cm to 
lxlO 12 cm, (values obtained from Ertan et al. 2007) the 
range of Keplerian frequencies (r 3 K = GMns/47t 2 ^-) that 
should appear in the continuum noise is then from 0.058 Hz 
to 2xl0~ 6 Hz. This range overlaps with the detected range 
of broad-band noise in our work which supports that we 
could be detecting noise from the fall-back accretion disks 
suggested to exist in these systems (Ertan et al. 2007) . How- 
ever, the variations of this noise does not resemble the case 
for X-ray binary sources (see Figure 5). In X-ray binaries 
outburst effects last for a few days to a week (where one 
detects a change in the broad-band noise level, as well). The 
persistent level of emission changes with accretion rate from 
high states to low states of X-ray binary sources while one 
also detects reallocation of broad-band noise from high fre- 
quencies into low frequencies. On the other hand, in our 
entire data set we did not detect any reallocation of noise, 
the noise we detected always appeared in the 0.05-0.005 Hz 
range and our timescale for increased noise or high level of 
broad-band noise lasted for over a year or two which is not 
a detected phenomenon in X-ray binaries. 

As discussed in the above paragraphs we favor that the 
changes in the band limited low frequency noise in mag- 
netars are related most likely to the burst characteristics. 
Glitches seem to be the least effective of all in the broad- 
band noise changes from these sources. As discussed in the 
above paragraph, one origin for the detected broad-band 
noise in AXP and SGRs could be the suggested fall-back 
disks around these objects with the aforementioned caveats. 

The origin of this noise could also be partially related to 
the frequency fluctuations of the neutron stars as in isolated 
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Figure 3. The averaged PSD in time showing the broad-band 
noise of IE 1841-045 (top) and RXS J1708-40 (bottom). 
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Figure 4. The time averaged PSD of SGR 1806-20 showing the 
broad-band noise. 



pulsars and also to the variability that arises from the up- 
scattered radiation in their magnetospheres, particularly en- 
hanced in the presence of a magnetar corona (Beloborodov 
& Thompson 2005, Beloborodov & Thompson 2007). The 
plasma coronae in magnetars are a result of the starquakes 
and twisting of the magnetosphere that injects an electric 
current (i.e., plasma) into the stellar magnetosphere (Be- 
loborodov & Thompson 2005, Beloborodov & Thompson 
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2007). Emerging currents are maintained during the X-ray 
outburst and a dense thermalized plasma is present in the 
magnetosphere (lasting from 1-10 years). Broad-band noise 
due to processes such as photon delays between soft and 
hard X-ray photons are well studied in X-ray binary sources, 
however such timescales are much shorter like a few millisec- 
onds as a result of Compton-up scattering (in the presence of 
Thompson optical depth) from disk coronae (e.g., Vaughan 
et al. 1994, Schulz & Wijers 1993) which yields a contin- 
uum noise in very high frequencies in comparison with the 
frequency range we detect here (within 256 sec). However, 
the processes in the magnetosphere of AXPs and SGRs are 
more complicated. Lyutikov & Gavriil (2006) shows that 
magnetospheric plasma can distort the thermal X-ray emis- 
sion emerging from the star surface by repeated scatterings 
at the cyclotron resonance (e.g., resonant cyclotron scatter- 
ing). 

Recently, Rea et al. (2008) showed that an X-ray spec- 
tral model that involves resonant cyclotron scattering (to- 
gether with a blackbody model) successfully fits the soft X- 
ray emission (1-10 keV data) of ten magnetars. If the domi- 
nant scattering process is the resonant scattering rather than 
Thompson scattering (as shown by Rea et al. 2008), then 
since the Thompson optical depth to scattering is different 
than resonant scattering optical depth, the escape times of 
photons from the magnetosphere or an existing magnetar 
corona will be different in comparison with a simple disk 
corona where Thompson scattering optical depth dominates 
(as photons are Compton up-scattered in the corona). Once 
escape times differ, so will the time delays between soft and 
hard photons will be different in a resonant scattering dom- 
inated plasma. Rea et al. (2008) calculates the ratio of op- 
tical depths as r res /rT = 10 5 which strongly suggests that 
given similar geometries the high level of resonant optical 
depth to scattering should cause more delays in the emerg- 
ing photons. If one assumes that scattering optical depth is 
directly proportional to number of scatterings, a millisecond 
delay then should become about 100 sec delay, that could 
yield noise around 0.01 Hz. We suggest that especially for IE 
2259+586 this scattering can produce or contribute to the 
band-limited red noise we detect in the 0.005-0.05 Hz band 
and its changes. This is, also, supported by the commensu- 
rability of 2 years rise in noise we detect in two AXPs and 
the magnetar coronal time scale of 1-10 years (time scale the 
magnetar corona exists). Moreover, the increase in the pho- 
ton indices of power-law tails in the X-ray spectra: For IE 
2259+586, V increases from -4.3 to -2.8 after the outburst 
(Woods et al. 2004, Kuiper et al. 2006); for IE 1048.1-5937, 
r increases from -3.4 to -3.0 slowly after the reported burst 
at 53185 MJD and this increases to -2.3 after the flux out- 
burst in 2007 March (Tiengo et al. 2005; Campana & Israel 
2007, Rea et al. 2007, Tarn et al. 2008). 

Characteristic of the activity is important for inducing 
such a process, e.g. a long burst or an outburst like char- 
acteristics are necessary to produce such effects of highly 
increased broad-band noise levels. Glitches or flares show 
no effect or slight variations in noise levels. This supported 
by nearly stable band limited noise levels of 1RXS J1708-40 
and IE 1841-045 which show only glitches and no bursts 0. 



2 http://www. physics. mcgill.ca/^ pulsar/magnetar/main.html 



Moreover, isolated pulsars that show glitches do not show 
such long term significant rises in red noise levels in the low 
frequency band. 

Other models explaining the magnetar timing variabil- 
ities are ruled out mostly because of their time scales. For 
example, seismic vibrations that are used to explain QPOs 
are very rapid and connected to giant bursts, and thus is not 
a satisfactory explanation for broad-band noise. One mod- 
ifying principle for the seismic models is the propagation 
of Alfven waves towards the core. Levin (2006) noted that 
these oscillations are a continuum, but they decay rapidly 
thus effectively they can not correspond to our persistent 
rms noise observations. 



CONCLUSIONS 

We observed broad-band noise (BLN) in the four AXPs (IE 
1048.1-5937, 1RXS 1708-40, IE 1841-045 and IE 2259+586) 
and an SGR (SGR 1806-20) mainly over the 0.005-0.05 Hz 
frequency band. The integrated rms values were in a range 
2.5-15% for four sources and 0% rms for IE 2259+586 before 
bursts. After bursts (our data do not include burst /outburst 
epochs) 20-70% integrated rms noise was observed in two 
cases; IE 2259+586 and IE 1048.1-5937, respectively. The 
broad-band noise characteristics do not resemble either the 
accreting X-ray binaries or the isolated pulsars completely. 

In general, we found that the four sources indicate a 
persistent level of integrated rms in time over a long base 
line of 3.5-5.5 years which may be affected by noise from fre- 
quency fluctuations, and most likely as a result of variability 
due to resonant cyclotron scattering in a magnetar coro- 
nae and/or variability from existing fall-back disks around 
these objects. We also detected large scale changes of the 
broad-band noise particularly due to long burst/outbursts. 
SGR-like bursts did not produce these type of variations in 
broad-band red noise, however we believe that this is lim- 
ited with the sensitivity of RXTE (also its background) and 
statistical quality of RXTE light curves. We calculate that 
IE 2259+586 displays high level of broad-band noise lasting 
for a timescale of 2 years following its outburst and glitch. 
In addition, we observed significant rise of the rms broad- 
band noise level in IE 1048.1-5937 for about 1.95 years, after 
the burst epoch at 53185 MJD. We believe that this noise 
can relate to the existence of remnant disks suggested to ex- 
ist around these objects, but the timescale of changes after 
burst/outbursts are inconsistent with the observations in X- 
ray binaries. We stress that the changes in the broad-band 
noise levels are related to resonant cyclotron up-scattered 
radiation in AXP and SGR spectra. Correlation with long 
bursts and outbursts indicates an origin in the variability as- 
sociated with the magnetar coronae. The typical timescales 
of existence of magnetar coronae of 1-10 years is consistent 
with the long term variation of BLN detected in our work. 
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